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Abstract 

In this paper, we present an improvement of a method for computing scattering amphtudes that 
include external (polarized) fermions with the following features: the formulas are quite general 
and work for different kinematic configurations and different external masses, they are explicitly 
covariant, they do not depend on a specific representation of the Dirac matrices and they have a 
meaningful limit when the masses tend to 0. The results presented make use of some well known 
formulas describing the density matrices in terms of projection operators within a more general 
formalism. Since our formulas intend to be as general as possible, we also take into account the 
possibility of a transverse polarization for massless fermions. 
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The notion of spin is somewhat difficult to express in a general, covariant formalism 
which would have also a meaningful limit for m — 0, since the Poincare little groups are 
different for > and = 0. Therefore, for the computation of Feynman amplitudes, 
we cannot obtain general formulas that are exactly continuous in the limit m ^ 0. 
Furthermore, when fermions are involved, one has to deal with spinors and Dirac matrices, 
and no particular representation of these matrices has some physical meaning. But for 
practical calculations involving spinors, a lot of people still make reference to a specific 
representation of these matrices, which leads to analytical results that are not explicitly 
covariant. Therefore, for the explicit calculation of some elements of the S matrix which 
involve fermions, there could be an improvement if one can express them within a very 
general, covariant, representation independent formula, where the spin degrees of freedom 
clearly appear, and where the limit m — > is meaningful and rather straightforward. 
There already exist in the literature some very useful formulas (see ||, |^), and the 
purpose of this paper is to generalize them, keeping in mind that we want to respect the 
criteria described above. For instance, the definition of the "conjugate" spinor ip by ?/^'^7° 
is purely conventional , assuming this way that we take a unitary representation of the 
Dirac matrices. This actually is not imposed by a physical principle but is rather a way 
to normalize the lagrangian. These matrices are even not supposed to be hermitian or 
anti-hermitian, and most of the calculations in this paper are completely representation 
independent, except for the sake of illustration. 

It is not claimed however that the formulas presented in this paper lead to faster 
algorithms for the calculation of huge Feynman amplitudes, and we must recall the reader 
that for this purpose, there are some quite fast calculation techniques ^ using spinor 
inner products. However, these algorithms are very fast only for massless fermions and 
introducing the case of massive fermions requires much more complex calculations. It is 
possible that in some cases, the formulas presented in this paper could be competitive 
for writing Monte-Carlo programs. Our purpose here is mostly to give simple tools for 
calculating analytically some amplitudes that are not very large, in such a way that 
one can possibly see the main physical features of the amplitude just by looking at its 
expression. For spinor inner products, one makes use of polar coordinates which need to 
set several geometrical conventions, yielding calculations that are not explicitly Lorentz 
covariant. This is why, although these tools exist and have proved to be efficient, we have 
found also useful to present how to calculate some amplitudes in a way which respect the 
symmetries (physical or not) of the problem and which includes as few conventions as 
possible. 

The paper is organized as follows. In the ffist section, we review some generalities 
about the Dirac equation, focusing especially on the subtleties that are very scarcely 
found in elementary textbooks, and which we found useful to gather here. Then we show 
on a simple example of calculation involving neutrinos, how one can in general express an 
amplitude in a representation independent and explicitly covariant way. The next part is 
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devoted to the mathematical derivation of our general formulas. For this purpose, we will 
redemonstrate the well-known formula giving the density matrix of a pure spinor state 
uu = (J) + m)(l + ^■y^)/2, using some representation independent calculations. We then 
show how to compute Feynman amplitudes or its square using these projection operator, 
focusing especially on the possible singularities appearing in the phase space which must 
be taken under consideration for an implementation in a Monte-Carlo program. The ad- 
vantage of the demonstration we used is that it shows the uniqueness of the decomposition 
of the density matrix, it exhibits clearly the spin degrees of freedom, and especially the 
transverse degrees of freedom naturally emerge for massless fermions. 



1 Generalities about the Dirac Equation 
1.1 From Klein-Gordon to Dirac 

In this section, we will review some basic things about the Dirac equation, and we will 
also focus on some points that are present in the literature, but unfortunately in very few 
papers. 

First, one must recall that the Dirac equation is obtained by looking for a factorization 
of the Klein Gordon equation: 

[d^g^''d, + m^]^ = (1.1) 
provided we choose the (H ) metric convention, and which can be factored into: 

(i^ + m)(z^-m)^ (1.2) 

where ^ = 9^7^, the 7^ are the Dirac matrices (4x4 complex matrices which must 
obey the anticommutation rules {7^^, 7''} = 2g^^). Our choice of a specific sign convention 
for the metric is of importance here. If we choose the ( — h + +) signature, we must then 
replace m by im or, to keep the mass term real, 7^ by Z7^, which is not algebraically 
equivalent to 7^^. This could lead to some confusions when one reads the large literature 
about the differences in the structure of the two Clifford algebras corresponding to the 
two possible sign convention for the metric, (1, 3) or (3, 1) (see and also some more 
mathematical references P, |10|, [ill]). However, the resulting physics fortunately does 



not depend on this metric convention, and to show this, we shall include the two possible 
sign conventions by rewriting the Klein- Gordon equation in the following way: 



[r]d^g^''d, + m^]^> = Q (1.3) 
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where rj = +1 if the signature is (1, 3) and rj = —1 ii the signature is (3, 1). Thus the 
factorization now reads: 

{iuj(^ + m){iuj(^ - m)-^ (1.4) 

where u is one square root of rj, i.e. up to a sign. Then, we can also choose between 
the two operators + m) and — m) to define the Dirac equation which gives 
another sign ambiguity. Therefore, to reabsorb these ambiguities, we define the matrices 
7^^ = 0^7^. We will consequently denote v = v^'^^ where is any 4-vector. The new 
anticommutation rules are therefore {7'^,7''} = 27] g^^ . That is to say, the 7^^ matrices are 
in a representation of the C(l,3) Clifford algebra whatever the sign convention for the 
metric is. The study of C(3, 1), with a completely different complex structure becomes 
therefore irrelevant on the physical point of view. We can now set the Dirac equation 
to be {id — m)\E' = 0, and since we have just shown that the only relevant signature is 
(1,3), we will assume throughout this paper that this metric is chosen (or equivalently 
that 7/^ = ^^'). 

Since the 7^^ matrices (or 7^^) do not transform like a four-vector, there must be a 
specific transformation for ^! to ensure the Lorentz covariance of this equationEH: 

*(x) 5a^(A(-^)x) = exp ^{{^^-^^Y,x'') (1.5) 

with a^" = f[7'^,7'^] (or equivalently f[7'^,7'^]), and Uf^^i, is an antisymmetric tensor, 
depending on the Lorentz transformation one performsE^. Thanks to this transformation, 
the Dirac equation {{id — m)\E' = 0) is Lorentz invariant. 

Now one can question how arbitrary are the Dirac matrices. First, it is easy to check 
that if we have a set of 4 Dirac matrices 7^^, and S is an invertible 4x4 matrix, then 
the new set 7"^ = S'y^S'^ still obey the anticommutation rules, and can be used as well 
as the first set of matrices, provided we operate the transformation \E'' = 5'\E' on the 
wavefunction. The very interesting property is that the converse is true, and we have 
summarized in the appendix the old proof given in Q. This property will be used in the 
following section, where we will define the notion of charge conjugation. For this purpose, 

(2) \Ye recall the reader that by Lorentz group is understood the subgroup of 0(1, 3) which is connected 
to the identity. It is also called the orthochronous Lorentz group and denoted L^^. Spinors are in a 
representation of its covering group SL{2,C) (see |l^): but are not in general in a representation 
of the full Lorentz group, covered by what one call the Pm(l,3) group (and Pm(3, 1) for 0(3, 1)). 
The fact that one consider only comes from the possibility of parity or time reversal violations 
(see m, ^ m, |9) and related experiments [|^ ||), but some authors argue that, for systems 
that are symmetric under the full Lorentz group, one can in principle see experimentally some differences 
between Pm(l,3) and Pm(3, 1) 

(3) 5*^ obeys S'[~^^7^S'a = A'^^j''. Using the expansion e^^Be-^ = B + [A, B] + [A, B]] + . . . and 
[cr^.^'^^p] = 2i{g''P-fi' - gf^P-f"), one gets this relation: A''^ = exp[-2(tj)]^^. 
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we will also need the general solution of the Dirac equation, expressed in momentum 
space: 



(27r)32/ 



{a{p)u{p)e-'P-'' + ¥{p)v{p)e'P-'') 



;i.6) 



where we have included the creation and annihilation operators {a{p) annihilates a 
fermion with momentum p and 6^ creates an antifermion) . u and v are the associated 
spinors that are respectively solutions of (jf^ — m)u{p) = and {"^ + m)v{p) = Q^. 



1.2 Spinors and Charge Conjugation 

The Dirac equation being given, we must now construct the corresponding lagrangian. 
For this purpose, we have to face the definition of adjoint spinors. Indeed, when one 
calculate a specific Feynman amplitude, one commonly use the definition ip = ?/'^7° and 
the property 7'^^ = 'y'^■y^^■y^. One must be careful to the fact that the use of 7° here is 
only conventional, and has very little to do with physics. 

The main point is to see that whenever {7^,/i G |0..3]} obey the anticommutation 
rules, then {j^^ fi G [0..3]} verify also the same anticommutation relations. We have 
then demonstrated (see in appendix) that these two basis are related by an interior auto- 
morphism, i.e. there exists an invertible matrix such that: 



V/i G I0..3] 7^*^ = H+rH+' (1.7) 

and also it implies (using the Schur lemma) that verify the property hI_ = h^H^ 
with h_|_ a complex number of modulus one. Then the adjoint spinor is defined a.s ip = 
ip'^H^. The Dirac lagrangian can then be constructed as follows: 



C = ^{i^-m)ij (1.8) 

For the sake of unitarity of the scattering matrix, we need the classical lagrangian to 
be real, which implies h_|_ = 1, that is to say if+ is hermitian^. For the class of unitary 
representations, we can take = 7*^. However, unitary representations have no more 
physical relevance than any other representation, and therefore we shall not suppose in 

^■^^We can remark that we could have chosen to write the Dirac equation with a +m instead of —to, 
which is equivalent to make the transformation ip j^ip on the spinor, and it is in fact just a change 
in the representation of the gamma matrices (because it is equivalent to make the transformation 7^ — s- 

(5)When one changes the representation of the Dirac matrices through 7^ — Sj'^S ^, changes into 
= XS''H-^-S like an hcrmitian form (its signature is (2,2)). 
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the following that iJ+ is equal to 7''. Yet, since if+ is defined up to a complex constant, 
imposing the unitarity of if+ is a good way to normalize the lagrangian. We will rather 
impose a normalization on the density matrices at the end of section [4 .31 . Although it 
leads to the same constraint at the end, it allows us to keep undetermined. In other 
words, we will normalize the density matrices instead of 

Now, the notion of charge conjugation needs the introduction of the gauge field and 
we get (i^ — ej^ — ■m)tp = for the equation of motion of the interacting fermion field. 
Taking now the complex conjugate we obtain the charge conjugated equation (z^ + 
— m)ijj'^ = 0, provided we take the following definition for ip'^: 



= Cl^^* (1.9) 
V^Gl0..3] 7; = -C-7,CZ' (1.10) 

(Note that when one conjugates ip, one shall not forget to take the hermitian conjugate 
of the creation and annihilation operators). This definition can be seen at the level of free 
spinors in the momentum representation: if we have — m)u{p) = 0, then after complex 
conjugation one gets + m)Cz^u* (p) = 0. In other words, the charge conjugated spinor 
of a so-called u spinor becomes a v spinor, which is the intuitive definition of charge 
conjugation (by computing tp'^ one can make the following identifications: u'^ = Cz^v* 
and v'^ = Cz^u*). The existence of C_ is actually given by the theorem we have proved 
in appendix relating sets of matrices that obey the anticommutation relations, which is 
the case of — 7*!^. C_ must satisfy C* C_ = c_/ where c_ is a real constant^. 

Having defined charge conjugation, we can take the opportunity to define a Majorana 
spinor. A Majorana spinor is its own charge conjugated particle, which necessarily con- 
strains this kind of particle to be neutral, and it must therefore obey the relation ip'^ = aip 
{a being a possible phase). It implies the following relations between creation-annihilation 
operators and u and v spinors: 



^^•'In the so-called Majorana representation of the Dirac matrices, C_ is just a multiple of the identity 
matrix, and if we look at another representation 7^'^ = Sjf^S^^ then the matrix C_ transforms with the 
formula C'_ — aS*C-S~^ , where a can be any complex number different from 0. 

(^)ln fact, one can scale the matrix C_ such that c_ ~ ±1. The sign of c_ depends on the convention 
we choose for the metric, but not on the representation of the Dirac matrices. Actually, one has c_ = +1 

for the convention (H ) (see |^) and c_ = —1 if one uses the ( — h + +) convention, like in the 

book of Jauch and Rohrlich Q]. In this latter case, charge conjugation becomes an anti- involution, and 
we could wrongly state that Majorana spinors (spinors which are identical to their charge- conjugated 
version) do not exist in this case, and the physics would depend on the convention for the metric. The 
point is in this case that the rule given to define charge conjugation would change, and we should have 
with this convention ijj'^ = CZ^4' where 7* ~ C+7pC^^ and C^C+ = c+J (c^ > 0). One can obviously 
relate C+ and C_ by C+ ^C--^^. 



5 



b = a (1.11) 
CZ^v* = au (1.12) 
CZ^u* = av (1.13) 

And for the sake of consistency, one should have C^C^ = j^-^- This happens to be 
possible because the sign of c_ is independent of the representation of the Dirac matrices, 
and in the Majorana representation of the Dirac matriceJi^, C_ is proportional to the 
identity (and thus C_ can always be written in the form C_ = ^/cZS*S^^ where S is 
any 4x4 invertible matrix). In fact, eq. |1.12| and |1.13| don't tell us much more about 



the structure of the wave function of a Majorana fermion, and Majorana spinors are 
nothing but standard spinors. The main difference comes from the quantization of the 
field, and especially from eq. |1 . 1 1| which allows to construct specific lagrangians that 



may distinguish between Majorana and Dirac spinors. We won't give more details about 
Majorana fermions here and refer the reader to the literature on this topic. 

In a similar way, we could also have defined an invertible matrix T+ or T_ for the 
transposition of the Dirac matrices (which also obey the anticommutation rules), but the 
three operations are related and since T_ can be obtained by combining if+ and C_, we 
will not have to use it. However, it can be useful to get a relation between C_ and 
using the fact that hermiticity and complex conjugation are two commuting operations. 
For this purpose we write: 



\ = illr = -HIC^J,CZ\HI)-' (=-T_7,T_-^) (1.14) 

% = ii;y = -{czyH^^,H-'ci (i.is) 

And using the Schur lemma, we know that there exists a non-vanishing constant k 
such that: 



= K{CzyH+ (1.16) 
H-^Cl = kCz\HI)-^ (1.17) 

Then, from a density matrix like p = iptp (or more generally a rank 1 matrix like ipip'), 
we can compute its charge conjugated version i/j'^i/j'^, and using the relations written in 
this section we may obtain easily: 



(^^Note that the Majorana representation is only a kind of "eigenvector" for complex conjugation (up 
to a real inner automorphism), which is different from the notion of Majorana spinors. The existence 
of a Majorana representation (which depends only on the metric signature), only implies the possible 
existence of Majorana spinors. 
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Figure 1: Amplitudes involving some fermion pairs 




;i.i8) 
;i.i9) 



The last equation being obtained thanks to eq. |1.17| , and if we require the lagrangian to 
be invariant under charge conjugation we must set k = 1. We will see beyond (section I 



that charge conjugation, given by eq. |1.19| can be expressed as a condition on some Lorentz 
tensors which define completely the density matrix p. 



2 General fermionic amplitudes 

In this section, we sketch the way we intend to compute any scattering amplitude with 
at least one fermionic current. We shall first note that there are four different possible 
configurations corresponding to a priori four different expressions for the amplitudes (see 
fig. H, where the overlined letters correspond to the antiparticles) . 

In a condensed notation, one may note a general amplitude as: 

M = i)'Oi) = tr[0^^'] = tr[OS] (2.20) 
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With E = u{—cp)u{—c'p') where u{p) is a solution of — m)u = 0, and c, c' are the 
fermionic charge {i.e. —1 for a fermion, +1 for an antifermion), provided that p is the 
momentum corresponding to an incoming fermionic current, and p' is the momentum 
corresponding to an outgoing fermionic current (see fig |I]). Of course, we could have 
chosen to write an amplitude by reference to an antifermionic current, or an annihilation 
or a pair production, but the four type of amphtudes are related and it is just a matter of 
convention. Also, we will parameterize the generic amplitudes using u{—cp) for the two 
possible choices (fermions or antifermions) instead of using charge conjugation because 
we will see beyond that charge conjugation can be easily expressed in terms of changing 
the sign of the constant c or c' and also the sign of the spin degree of freedom which 
will be introduced in the next section. The general expression of the amplitude can be 
specialized to the four kinematic configurations of figure |l] in this way: 



= vip\)Ouipl) = tr[Ouipl)ui-p\)] (2.21) 

M^n,out = u{p'L)Ou{pl)=tr[Ou{pl)u{p'L)] (2.22) 

Mj;^^^, = v{pl)Ovip\)=tr[Ou{~pl)u{-pl)] (2.23) 

A^o«*,^ = u{p'L)Ov{pl)=tr[Ou{-pl)u{p°,)] (2.24) 

where p_ denotes the momentum of the fermion, is for the antifermion, and O is 
the interaction operator. Instead of decomposing spinors in their components, we will 
express the amplitudes without having to specify a specific representation for the Dirac 
matricesElll. We shall therefore try to express the rank 1 matrix S = u{—cp)u{—c'p') in 
the basis of Dirac matrices, and we can do the same for the operator O. Then the trace 



can be easily computed using simple algorithma^MI. If one prefer to compute directly the 
square of the amplitude, one has: 



|>i|2 = [u{-c'p')Ou{-cp)][u{-c'p')Oui-cp)]^ 

= [u{-c'p')Ou{-cp)] \u{-cp)H-^O^H+u{-c'p'] 



tr 



p'OpO 



(2.25) 



u{—dp')u{—c'p'] 



and O = H.^O^H^. If we decompose 



with p = u(—cp)u(—cp), p' 
O on the basis of the Dirac matrices we have: 

'^^■'Some variants of this formulation exist already in the literature, see ||, ^. Another calculation 
technique was developed by Hagiwara and Zeppenfeld decomposing spinors in the chiral represen- 
tation into their Weyl spinors. This method yields a quick algorithm for computation but here, we are 
looking for a formulation that is explicitly covariant and representation independent. 

Symbolic calculation programs can do this very easily nowadays. The fastest being probably FORM, 
and a rather convenient one is "M" iBSl . 
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(2.26) 



an 



= s:- sl^' + C + Kl' + S%*a'^-' (2.27) 

Then, in the following sections we will give the decomposition on the same basis for 
the density matrix p, which will give the final result for E and Ai. Before, we shall discuss 
a simple example to show the practical procedure we will use in the general case. 



3 A simple example 

We will thus consider in this section the scattering of a neutrino, for instance in an 
electroweak process. The neutrino spinors can be easily expressed in terms of Weyl spinors 
if we take the chiral representation for the Dirac matrices. We get: 

-e-^/2sin(^/2) \ ^ _( cos(^/2) \ 

e*^/2cos(^/2) j ' ~ 1^ e^^/2sin(^/2) ) ^'^'^^^ 

and using the notations P[in,out) = -£'(m,oiii)(l, ''^(in,oMt))) the S matrix can be computed 
with a little bit of algebra: 




out J 



Q 




EinEoutl 





VL 

1-^5 



(3.29) 
(3.30) 



One may then define the vector such that f° = (i7ii + fi22)/2, = — (i7i2 + ^^2i)/2, 
v"^ = —i{Qi2 — fi2i)/2, = —(fill — ^22)/2. When the two fermions are not back-to-back, 

doesn't vanish and one can impose the condition that is real, multiplying by the 
suitable phase. One finally obtains the following form for f^: 



EinEgut 



V2 



^in ' ^out 



^out 



iP'in + ^out + i^in A ^out) 



\ 



J 



(3.31) 



One can see that is orthogonal to both momenta of the fermions, and moreover, 
that the space-like part of this vector represents an elliptic polarization associated to the 



(11) 



By we denote v*^^^ . The Dirac matrices are not conjugated in this notation. 
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3- vector nj„ — Uout- We have therefore combined the two spin (1/2) particles to obtain a 
spin-1 current. 

But we will show that the expression of w'^ can be expressed in a covariant way through 
the formula: 

_ (Pout ■ k)P^n ~^ (Pin ' k)Pout ~ (Pout ' Pin)^^ ~ '^^'^^'"^^i^Pin. pPoMt ct 
\f^\ll{Pout ■ k){pin ■ k) - {pout ■ Pin)k'^ 



Where the four-vector k^^ serves as a reference. Expression p.31 is actually obtained 



by setting k^ = (1,0), but one can note that is invariant under the transformation 
k ^ k + Xpin + fJ'Pout, which also implies that k must not lie in the plane generated by 
Pin and Pout- This will be rather obvious with the following. Note also that a rescaling of 
k does not affect v. Therefore, there is only one degree of freedom of k that may influence 
V, through a phase multiplication. 

To obtain this expression, let's consider the density matrix for both the incoming and 
outgoing momenta (we will derive its expression in the following section): 

1 + 7^ 

uu = -f) — ^ — (3.33) 
then one has for a given vector k^: 

1 + 7^ 1 + 7^ 

UinUin^UoutUout = (Uin^Uout)^ = An^^out ^ = ^ (3.34) 

with: 



V'' = {Pout ■ k)p^n + {pin ■ k)p^out - {Pout ■ Pin)k^ " i^'"''"' k^Pin pPout a (3.35) 

If k lies in the plane generated by Pin and Pout, v vanishes and the latter relations are 
useless. So let's take k outside this pla ne. Then we obtain the normalization factor in 
eq. 3.34] by multiplying once again by and taking the trace, we get: 



,1 + 7^ 



tr 



2{v ■ k*) 



2 

2 (23?((pi„ • k){pout ■ k*)) - {pout ■ Pin){k ■ k*) 
+ie^'"'''k^klpinpPouta) (3.36) 



(12)^* fc*'^7p (see previous section). 
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e''^V2^2^{ipi^ ■ k){pout -k*))- {pout ■ P^n){k ■ k*) + ie'^'^P'^kf^k 

uPin pPout a 

(3.37) 

In this example, we have been obhged to introduce another 4-vector which has a priori 
no physical meaning (its role is only to set a reference for the phase of the amplitude), and 
we can be led into trouble if this 4-vector is not suitably chosen. This problem would not 
occur if we had computed directly the square of the amphtude, for which the introduction 
of this arbitrary 4-vector wouldn't be necessary. This also enlighten the impossibility to 
find a formalism for the amplitude which would be continuous in the massless limit, since 
this kind of singularity does not exist if one of the fermions is massive and k is chosen 
on the light cone. We will discuss this point further at the end of the paper, and we now 
turn to the general case of an amplitude that includes a fermionic current, represented by 
what we have called the S matrix. 



4 The general case 

4.1 General properties of the S matrix 

In this section, we will specify the properties of the S matrix. Instead of expressing the 
spinors in a specific representation of the the Dirac matrices, we will express the general 
equations that determine S without having to take any specific representation. These 
equations are therefore more physically meaningful. This kind of formalism have already 
been studied a very long time ago by Pauli and others. It was called the bispinor algebra, 
and their intrinsic relations were studied in the context of Fierz identities (see p4|). We 
just give here another demonstration, which does not use Fierz identities which contain 
far too much information for our purpose. The main point of our proof being in fact a 
representation independent characterization of a rank-1 4x4 matrix. To proceed in the 
same way as in the previous example, we shall deal first with the density matrices p, for 
which we will use three constraints. The first is the equation of motion {'p — m)p = = p{'p 
— m). Then we will use a conjugation constraint p H^^ = p, and finally a 
criterion to say that p is a rank 1 matrix (using a characterization described in appendix). 
In principle, after these three steps, there should only remain a normalization condition 
to determine. So, let us first write the consequences of the equations of motion on p, 
provided we have decomposed p on a typical basis of Dirac matrices, that is: 

p = s + s^^ + ^ + 4^^ + S^.a'^" (4.38) 
where S^u is an antisymmetric tensor, a^'^ = |[7'^,7'^], and one has the relations: 
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s = ltr[p] s = itr[p-f^] 
Then the equations of motion can be also projected on the same basis and it gives 



{j>-m)p = 0^ (4.39) 

p ■ V = ms (4.40) 

p ■ a = ms (4-41) 

spf" -2iS^''p^ = mv'' (4.42) 

sp^" -e'"^^"'Sai3Pu = ma^ (4.43) 

\E^-''-apP^-\\p,vf^ = mS'^" (AAA) 



and: 



p{j()-m) = ^ (4.45) 

p ■ V = ms (4.46) 

p ■ a = —ms (4.47) 

spf" + 2iS'"'p^ = mvf" (4.48) 

-sp^" -e^^f'^S^f^p^ = ma^ (4.49) 

\e^^("^apP,\\\p,vY' = mS'^'^ (4.50) 



These two sets of equations can be rearranged into a simpler one: 



s 

[p,vr 

sp^ 
p ■ a 



(13) "^il pf^v" 





O^v^" = Ap^ 
mv^ => s = Am 





-ma'^ 



2mS 



(4.51) 
(4.52) 
(4.53) 
(4.54) 
(4.55) 
(4.56) 
(4.57) 
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In the massive case, the two last equations are indeed equivalent. We can already see 
that the structure of the density matrix is quite constrained. Since s = 0, the conjugation 
constraint simply says that all the remaining coefficients {s, v^, and Sfj_y) must be real 
(see eq. |2.27| ). We shall now use the characterization of p as a rank 1 matrix, which is 
the last constraint on p. In order to find the necessary and sufficient conditions for p 
to be a rank 1 matrix, we will apply the corollary of the second theorem demonstrated 
in appendix B, and in the ffist place we shall not make reference to the constraints we 
have already obtained with the equations of motion. We will add them at the end. The 
corollary tells us that we must have tr{pQpQ') = tr{pQ)tr{pQ') for every matrix Q and 
Q'. We will therefore choose for Q and Q' the different kind of Dirac matrices, which 



leads tdlM 



-k For Q = I we obtain respectively for Q' = I, 7^, 7^, 7^7^, a'^'^: 

2Sa(3S'"' = 3s^ - - {v^ - a^) (4.58) 

2ss = -ze>^-P-S^,S,^ (4.59) 

sv^" = e^''^^'axSaf3 (4.60) 

sa^ = e^'^'^^'vxSaf^ (4.61) 

25^^^^ = -e''^^"'{vaap + isSap) (4.62) 

Note that if s 7^ 0, the lemma demonstrated in the appendix tells us that these 5 
conditions [4.58| to |4.62| are sufficient to ensure that E is a rank 1 matrix. However, to 



treat also the case where s = 0, we must set Q to the other kind of "Dirac matrices". 
-k For Q = 7^ we obtain {Q' = 7^, 7^", 7^7^, a^"^): 

2^„^^"^ = 3f - + {v^ - a^) (4.63) 

svf" = 2za„5"'^ (4.64) 

sa^" = 2waS''f' (4.65) 

+i(S'^„S'p.£'^-^'^ - S\S,,E^^n (4.66) 
-k For Q = 7" and = 7^^, 7'^7^, we get: 



= ^"/3(_a2 - + - + 2^p,S''") - 2t;V + 2a°a^ + 8^";,5^^ (4.67) 

= c,"/3(a.t;) + (t;"a/5-t;V)-se"^^"S'p, + 2is^"'^ (4.68) 
= isia'^g'^'' - a!" g""^) + soa^^"'''' 

+2t;A(^^V - 5^"^°^) + 2(-^;°^'^" + - v"" S""^) (4.69) 
(i4)^g ^gg jjgj-e the convention £0123 = 1- 
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For Q = 7^7^ and Q' = 7^7^, cr^", we obtain: 








^"'^(-a^ -v^-s'^ + s'^- 25^,5"") + 2v''v'^ - 2a"a^ - SS^S^f^ (4.70) 
+2aA(^"^5^^ - ^"^S^'^) + 2(a"5^'^ - a^S"'^ + a''5"'^) (4.71) 



* For Q = a^^ and Q' = a^"', we obtain: 



^gf^Pg'^- - ^^"^''^(s^ + s2 + _ ^2 ^ 25„^5"^) 

+2 (^'^^(a'^a'" - t;''^;'") - ^^'^(a'^a^ - t;^^;'') - ^'^^(a'^a'" - v^v"") + ^'''"(aV - f^f^)) 



(4.72) 



In eq. 4.69 we can insert eq. 4.65 and it simplifies into: 



And similarly, inserting eq. [4.64| into eq. [4.71| we get: 



(4.73) 



(4.74) 



Realizing the sum of eq. [4.67| and eq. [4.70| we obtain = —v'^, which then leads to 



2Sai3S^^ = 2(s^ + s^) and f ^ = — using eq. [4.58| and eq. [4.63| . The latter relations 
can be inserted in eq. [4.67| and eq. |4.70| , the difference of which leads now to the simpler 
relation: 



Inserting this equation and the previous ones in eq. 4.72 we obtain: 



(4.75) 



.ss 



(4.76) 



4.2 Solutions of the rank 1 equations 



If we have the conditions s 7^ and — 7^ 0, the solutions of equations |4.58| to [4.62 
are given by: 
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a ■ V 



2 -2 
S — S 



2(s2 - S2) 



The dual of S'^'^ then reads: 



(4.77) 
(4.78) 

(4.79) 



S 



-e 



pupcr 



5, 



pa 



(4.80) 



However, the relation [4.77| must always be true thanks to equation [4.68| . The sum of 



eq. [4.67| and eq. |4.7CI| yields = — f ^ in any case, and the difference between eq. |4.63| and 



eq. |4.58| shows that the fundamental relation |4.78| must also be true in any case. Note 
also that in eq. [4.66| , the last term could be surprising, but using [4.79| , one finds that its 
contribution is 0. The remaining part of this equation is actually the "dual" equation of 
equation i.62[ 



Now what happens if s = ±s. We can rewrite equation [4.62| using the new notations 
S± = S ± S* and lj± = s ± s: 



VpaA 



(4.81) 



Then if loj^ or lo^ vanish, the corresponding self-dual tensor is left undetermined by 
this equation and a constraint appear between v and a. 



4.3 Solution with the whole set of constraints 

At this stage we can summarize all the constraints we have obtained. The equations of 
motion for the density matrix led us in the massive case to: 

p = A(m + 7^) + + ^e^^'^'^a^VfiOpv (4.82) 

with a • p = and the rank 1 characterization only adds the condition a} = —\^m? 
because of course, a lot of the previous equations are redundant. The last two equations 
show that the 4- vectors a/X ± p lie on the light cone. That is to say we can write 
a = X{p — ki) = \{k2 — p) with kf = = 0, ki ■ p = k2 ■ p = m?- Actually, only one of 
these two vectors is useful since, if ki is given, one can take 2p — ki for k2- Therefore, let 
us consider any non- vanishing 4- vector k lying on the light cone. We can therefore set: 
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a'^ = hX (^p" - ^k"^ (4.83) 

Where h = ±1. The scale of k is of no importance in this expression, whereas its 
direction becomes a good candidate for defining the direction of reference for the polar- 
ization of the fermion. We can also look at the light-cone limit of this expression. If the 
space-like projection of does not tend to the one of k'^, the limit is obviously the limit 
of hXp^, and the possible singularity occurs when p becomes almost collinear to k. In 
order to study this limit, we can set p^^ = k^ + e5^ with e — (we must have chosen k 
suitably here, and especially /c" > 0). The expression of then reads: 

a^" ( \ 

= -A;M 1 + e—- + eb^ (4.84) 
h\ \ b ■ k j 

And the limit for exists (a^ ~ —hXp). Rewriting the whole density matrix we get: 



= x(m + ^ + h(^^-^^^^'' + mh(^l-^^^''^ (4.86) 
= X{^ + m)(^l + h(^l-^^^'^ (4.87) 

We recover for the massive case the well known result (p = + m) with a-p = 
and = — 1), in a slightly different form. This last formula is already very often used in 
the literature, and one may question about the interest in redemonstrating it. The first 
interest is to use this formula (especially in the form of eq. [4.85| ) to compute scattering 



amplitudes or the modulus square of these amplitudes with external fermions, (as we will 
see in the following) in a way that respects the criteria we have defined in introduction. 
And also, it shows naturally how the spin degrees of freedom emerge and the possibility 
of some transverse degrees of freedom in the massless limit. In fact, when one looks at 



the limit m — > in eq. |4.85| , the last term vanishlil^, and if p ■ k doesn't tend to 0, the 



term containing 7^ becomes equivalent to 7^7^, whereas when p- k also tends to when m 
vanish (i.e. p^ —>■ Ck^), then the 7^ term is equivalent to The m — > limit allows 

us to recover the two possible density matrices used in the literature for massless fermions 
(p = ^{1 ± 7^)/2). However, in the massive case we had more freedom for this density 
matrix because the direction of k is free. In some sense, this freedom seems "frozen" in 



(i5)T]^is limit is obvious if p ■ k does not tend to 0. Otherwise, if p • /c — ^ when m — s- 0, one has 
p'' — > C/c''. Then one can see that the hmit is by taking for instance = {E = sjnfi + p^, 0, 0,p) and 
fc — (1, 0, 0, 1). Then e^^°^(pofc3 — Pa^o) = E — p = p ■ k and the hmit is clearly also in this case. 
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the limit m ^ 0, but it is not really the case, and we shall study in details the solutions 
of the constraint equations for m = exactly. We can actually deduce from eq. [4.56| that 
the generic solution for S^''' when m = is S*^*^ = A[p, £:]^'^/2 with p ■ e = (eq. [4.551 ). 



Therefore, S^"^ is not necessarily 0, contrary to what we have deduced from the limit 
m —>■ 0. The rank 1 conditions then add another constraint, but relax also one. Since 
= in our case, the normalization of is no more imposed (the coefficient A in factor 
of eq. 4.83| is no more valid) and from eq. 4.57 one has only = A'p^ with a priori X' ^ A. 



From eq. |4.75| we get a constraint on the normalization of e^: 



We can therefore write the density matrix in the massless limit in the form: 



p = j^(A + AV + 2^) (4.89) 
= (l + hVl + e^-f' + if) (4.90) 

The emergence of the "transverse" degrees of freedom in the massless case is not really 
a surprise. Their origin comes from the fact that the little group of a massless momentum 
(/SO (2)) is of the same dimension (3) as the one of a massive momentum (S'0(3)), and we 
should therefore obtain the same number of degrees of freedom. However, the transverse 
degrees of freedom of massless particles are not observed in experiments, and we are led 
to discard them for neutrinos (see |^ eq. 2.5.38). In the context of the study of some new 



physics where some massless fermions not yet observed could exist, there is a priori no 
reason to discard them and it is also a reason why we present the most general situation 
in this paper. Finally, for the global normalization of these density matrices, we choose 
A = 1/2 to conform to a common normalization used in the literature. 

4.4 Charge conjugation in this formulation 

Since we have obtained the decomposition of the density matrices on the basis of Dirac 
matrices, we shall now study how charge conjugation operates on these decompositions. 



From what we have seen in section 1.2 we can write 



(4.91) 

s*J - 5*7^ - f + ^*7^ - Sl^a^"" (4.92) 
sJ - ^ + ^7' - 5/..^'^" (4.93) 
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The last equation comes from the hermiticity of p and the fact that we have s = in 
any case. Applying this on the explicit expression of p (eq. [4.85|) , the charge conjugation 
acts rather simply on p\ 

p{p, k, h) p%p, k, h) = p{p' = -p, k' = k,h' = -h) (4.94) 
or also this way: 

p' = -p, k' = 2p -k, h' = h] (4.95) 

p-k J 

We can see that it is equivalent to change the sign of h, keeping the vector of reference 



constant, or to keep h the same and to enforce a symmetry operation on k. Also eq. [4.94 



has the advantage to apply in the massless case where we can also write from eq. 4.90 



p(p, e, h) p\p, e, h) = p{-p, e, -h) (4.96) 



4.5 Computation of a generic amplitude 

We have now all the elements necessary to compute a generic amplitude of the form 
M. = tr[0'S\. We therefore wish to proceed in a similar way as we did in the simple 
example of section |^. Thus we can write the equivalence: 

, uuku'u' , ^ 

uv! = _ _ = (4.97) 



where the 4-vector k serves as a reference for the phase of the amplitude and the 
equivalence means "up to a phase". It is important to note at this stage that we can 
be led into trouble if an amplitude contains several identical particles. In this case we 
should have the same reference for the phase of two diagrams where two identical fermions 
are permuted, because of the possible interferences. In this case, this method shouldn't 
work because we have arbitrarily changed the relative phase between the diagrams. In 
this case, a covariant expression can be obtained through the computation of the square 
of the amplitude, which often means huge analytical expressions. We shall not discuss 
further this problem. Coming back to our single fermionic current, we can generically 
write the S matrix in the form: 



tr[pu^Pu4] 
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where one can replace pu and pu' by their expression given in formula [4.85| for massive 
fermions or eq. [4.90| for massless fermions. Then, using eq |4.98| , the amplitude can be 
written: 

M = ^rEMti (4.99) 

'tr[p^^p^>^] 



and one has to compute two simpler traces, instead of one possibly huge trace if 
one wants to compute the square of the amplitude using eq. 2.25| . As for the choice of 



k^, we have seen in the simple example presented in section |^ that we may have some 
problems if this vector is not suitably chosen, especially in the massless case, because of 
a kinematic singularity appearing in the plane generated by the two momenta. And if 
one intends to implement this formula inside a Monte-Carlo program, where the events 
are generated randomly, it is better to avoid this kind of singularity to be sure that 
none of the momentum configurations will be close to the singularity. Therefore in the 
massless case, one may choose k'^ such that k"^ = 1, and in the massive case, the explicit 



computation of the normalization factor in eq. [4.99| shows that one may preferably choose 
k'^ on the light cone. Since the explicit calculation of these traces can be done easily by 
some symbolic calculation programs, we shall let the reader do them if he is interested 
in such calculations, because it is often important to take advantage of the particular 
situations to choose k^, or the spin axes properly. 

We may end this section by mentioning that we can also choose a simpler operator 
than Ijt to be inserted into the normalization factor when at least one of the fermions is 
massive. We can simply use the identity matrix instead of l/t and we get: 

M = 'l^Md (4.100) 

Which is simpler to compute, and shall not lead to kinematical singularities in most 
cases. Since the denominator can be of order mm', it will be numerically better to use it 
for heavy fermions. 



5 Conclusion 

In this paper, we have shown a quite general method for the calculation of Feynman 
amplitudes or its square with external fermions. The mass of the fermions can be of 
any value, even in the same fermionic line. The formulas given are also independent 
of the representation of the Dirac matrices, explicitly covariant and have a meaningful 
massless limit. The formalism is not exactly continuous in m = in the strict sense since 
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there are some degrees of freedom that appear in the case m = that compensate the 
transverse degrees that are "frozen" into the two hehcity modes when m ^ 0. However, 
these transverse degrees of freedom are not observed for neutrinos and they are therefore 
not taken under consideration for these particles. For the study of some New Physics, it 
can be interesting to keep them. The other trouble with massless fermions comes when 
one computes a scattering amplitude with a current composed of two of these massless 
fermions. There can be some singularity in phase space for the normalization condition if 
one takes a fixed momentum for the reference phase, which can lay in the plane of the two 
external momenta for some kinematic configurations. This singularity does not exist for 
the computation of the square of the amplitude, which fortunately in the massless case 
may lead to expressions of reasonable size, contrary to the massive case. 

To obtain the general expression of an amplitude, we started from well known for- 
mulas giving the density matrices that we have redemonstrated in the most general case. 
Using this kind of demonstration, the possibility of a transverse polarization for massless 
fermions naturally emerges, and we have also characterized within this formalism the no- 
tion of charge conjugation and Majorana spinors. We expect that these results can be 
useful to get simpler analytic expressions for short Feynman amplitudes, or not-so-short 
amplitudes computed using symbolic calculation programs. We also expect a simpler im- 
plementation of fermions in Monte-Carlo programs, thanks to the fact that the formalism 
presented in this paper requires very few conventions, which was also one of the goals of 
this work. 
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Appendix A: some properties of the Dirac matrices 

Suppose that we have one set of matrices that obey the anticommutation rules. From 
this set of 4 matrices, one can explicit a basis of A^4(C)Ei^, 7*" (r G [1..16]) in this way: 



7^ = ^o"o(r)^^ai(r)^^a2(r)^^a3(r) (5.101) 

with a^{r) = 0, 1. We may also sometimes denote n^(r) = ao{r)+ai{r)+a2{r)+a3{r). 
When n^{r) is even, the matrix is said to be in the even subalgebra of the Clifford algebra. 



^^^'>M4{C) denotes the set of 4 x 4 complex matrices throughout the paper. 
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This basis is chosen essentially for the proof of the theorem relating two different sets of 
Dirac matrices (see beyond), but most of the time one uses t he m ore convenient basis 
(J = 7(,=i), 75 = Z70^i^2^3 _ 7V, ^'^^ = i[7^71)'^• 

With these definitions one can prove that tr[7r] = if 7^ is not the identity matrix 
(r 7^ 1) and that 7^7^ = '/7p,q7/(p,q) , where rjp g = ±1, and more importantly, for a fixed q, 
p I{p,q) is a permutation of |1..16] (and similarly for a fixed p, q ^ I{Pi<l) is also a 
permutation). Also rjpq depends only on the anticommutation relations and the metric 
convention, and not on the specific set (representation) of Dirac matrices chosen. Since we 
have also shown in the beginning of the paper that the only relevant signature is (1, 3), we 
also assume this metric to be chosen. We have also tr[^i'~^j\ = Cij, where the coefficients 
Cij vanish if and only if i 7^ j. 



Demonstration of the relation between different representations 

Suppose now that we have two sets of Dirac matrices 7^^ and 7"^ that obey the anticom- 
mutation rules. As before, we can define a basis for the 4x4 complex matrices 7.^, and 
another similar basis 7'^ defined in the same way but using 7"^ instead of 7^^. We shall 
also need another 4x4 complex matrix, F upon which we do not impose any constraint 
at this time. Then we define the matrix S as: 



S= Y: I'rFlr^-'^ (5.102) 
r=1..16 

We will show that S must be or invertible, and then that in this latter case one has 
7"" = SY (it is sufficient to check this property for r = /i G [0..3]). 

The proof comes from the following relation: 

^7* = E irFlr^-'ht 
r=1..16 

= E i iit,s)Fliit,sl-'ht (^ r^I{t,s)) 

s=1..16 

= E i i{t,s)Fls^~^ht''^htrit,s (^ def. of -fi^t,s)) 

s=1..16 
s=1..16 

= E itl'sFiy''^ (^ def. o/7/(t,s)) 

s=1..16 

= 7'*^ (5.103) 

Sometimes we will have to consider 7(r=i) which is just the identity, different from the usual 7(^=1). 
We will therefore write the "r = ..." when necessary, in order to avoid some confusions about the meaning 
of the indices. 
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This is the very fundamental property we need to prove our latest assertion. Now if x 
is a four vector such that Sx = 0, then every is in the kernel of S. It means that Ker S 
is stable through the action of the whole algebra, that is to say every 4x4 complex matrix 
(A^4(C)). If we consider an irreducible representation of the Dirac matrices, it implies 
that Ker S must be an improper subspace, i.e. or the whole spinor space C^. Therefore, 
S is invertible or 0. Now, why is the four dimensional representation irreducible? Suppose 
that K = Ker S is not 0, and let L be a subspace of such that = K ® L. Let 
be an endomorphism of such that the image of one (non-zero) vector in K lies in L (it 
only needs L to be different from 0). Now, the key point is that M. can be decomposed 
on the 7r basis. Since K is stable by all the 7^, we get a contradiction and L must be 0, 

We have now demonstrated that two sets of Dirac matrices are related by an inner 
automorphism if we can find 7^ 0. And it is possible to find F such that S is not 
0. Otherwise, for all x,y E C^, posing F = yx\ one has Y^Av^l' ry)^^lr'^~^^ = O5 thns 
J2r{y'^l'ry)lr''~^^ = and taking the trace, one obtains y^'j'f^^^i^y = Ay'^y which cannot be 
if F 7^ 0. So one can find F such that 5* is non-zero, then S is invertible and the two 
representations are related through an inner automorphism. □ 



Unitary representations 

In this paragraph, we will show that one can construct a unitary representation of the 
Dirac matrices from any irreducible representation. 

So, if we denote H = J^rlllr, one can use some arguments similar to the ones used 
in eq. p.l03| (with also the fact that 77^^ = 1) to show that ■jj^H'jfj^ = H. Now, it is 



clear that H is hermitian definite and positive, which allows us to write it as a "square" , 
i.e. H = h'^h and h is invertiblel^i^. Thus we have 'jj^h'^hj^ = h^h and therefore the 
representation given by h'jf^h'^ is unitary. □ 



Appendix B: mathematical proof of some elementary 
theorems 

In this appendix we will demonstrate two interesting properties. The second one is used in 
the body of this paper. The first one is rather trivial if one uses a specific representation 
for the Dirac matrices. The aim here is just to illustrate that most of the important 
properties within the Dirac formalism can be shown without the need of any specific 
representation for the 7 matrices. 

^^^^ft, is also unique if it is taken hermitian positive definite. 
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Theorem 1: Let p be a non vanishing four momenta. We denote m — if > 
and m — i\/—p'^ if p^ < 0. Let c he a constant which can he set to ±1 (if the convention 

for the metric is )), or ±i (if the convention for the metric is ( — \- + +), case 

which can he avoided if we replace j) hy p, as shown in the heginning of the paper). Then 
Ker(J) — cml) = Im{j[) + cml) . 

Proof: Let M4 be the four dimensional Minkowski space, and a four dimensional 

complex vector space, in which we define spinors. 

Since = c^m^J, it is clear that Im(J) + cml) C Ker(j) — cml). In the case m 7^ 0, if 
u G Ker(J) — cml), then -u = (j^ + cml) {u/ {2cm)), and the theorem is proved. In fact, 
the eigenvalues of are of order two, because has at most 2 non-vanishing eigenvalues, 
and its trace must be zero. Thus ;^ — cml are two rank two matrices, and we recover the 
elementary result that the space of states for a spin 1/2 particle is of dimension 2. This 
is also true when m = as we shall see. 

Now suppose that m — 0. We know that dim Ker^f) + dim Im^f) — 4, and we have 
seen that Im{f) C Ker(J)), and since p^ is a non vanishing vector, j) ^ ^ and one has: 
1 < dim lm{j)) < dim Ker^p). Thus we have two possibilities: 

• dim Im^f) = 2 = dim Ker{p). Then the theorem is proved. 

• dim Im(J)) — 1 and dim Ker(J)) — 3. 

We will therefore show that the second case is not possible. Let Ui, i e {0, 1,2,3} 
be a basis of C4, such that {ui,U2,Us} is a basis of Ker(J)) and Ui — '^Uq. The vector 
jf is of the form p(l,n) (with — 1). Let ft' and n" two vectors such that ft', ft" 
and n form an orthonormal basis of . Then if q' = (0,n') and q" = {0,n"), one has 
W ■ p) = W ■ p) = {q' ■ q") = 0, and ^'^ = g'^ = ±1 = j^"^. Therefore, the eigenvalues of 
^' and ^" are in the set {i, —i} or {1, —1} depending on the sign convention for the metric. 

One has Ui = "^{uq), then ^'ui = pi—^'u^), but if dim Im{p) = 1, we must have 
^'ui = X'ui with A' = ±i (or ±1), and similarly for ^"ui. Now ^'{^"ui) = \'\"ui = 
—^"{^'ui) = —X'X"ui. This would imply \'\"ui = 0, which is impossible. Thus the only 
possibility is dim Im{;p) = 2 = dim Ker{p), which impUes Im{][)) = Ker{p), and the 
theorem is proved also for m — 0. O 

Theorem 2: Let M he a n x n complex matrix, then one has: 

rank{M) = 1^ M ^0 and^Q e Mn{C) MQM = tr{MQ)M (5.104) 

Lemma: Let M he a n x n complex matrix, such that tr{M) ^ 0. Then : 

rank{M) = 1 ^ = tr{M)M (5.105) 
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Proof: 



If M is a rank 1 matrix, then one can write M — xy^ where x and y are two complex 
vectors. Since we have tr{M) —y^-x, one has = tr{M)M. 

Conversely, since tr{M) 7^ then X A (X - tr{M)) = 1, and the condition M{M - 
tr{M)) = implies that C" = Ker{M)®Ker{M -tr{M)I) {n = 4), thanks to the kernel 
decomposition theorem. Thus, using a basis compatible with this decomposition, M can 
be written in the form: 



/ 

tr{M) 



V 



tr{M) J 



(5.106) 



Where P is an invertible matrix. It leads to tr{M) — dim Ker{M — tr{M)I) x tr{M), 
thus dim Ker{M - tr{M)I) = 1 ^ rank{M) = 1. □ 



Proof of the theorem: 

If rank{M) — 1, then VQ, rank{MQ) — Q or 1, and applying the same reasoning as 
in the lemma, we easily conclude that MQM — tr{MQ)M. 

Conversely, if VQ, tr{MQ) = then M = 0. Therefore, in our case, we can find Q 
such that tr{MQ) 7^ 0. And since Q tr{MQ) is a continuous function, and the group 
of invertible matrices is dense in Ain{C), we can choose Q invertible. Then, multiplying 
on the right side by the matrix Q we get (MQ)^ = tr{MQ)MQ and the lemma tells us 
then that MQ is a rank 1 matrix, and thus M = {MQ)Q~^ is also a rank 1 matrix. □ 

Corollary: Let M be anxn complex matrix (M ^ 0). We have the following equivalence: 
rank(M) = 1 <^ V(g, Q') tr(MQMQ') = tr{MQ)tr{MQ') (5.107) 

Proof: 

If one knows that M is a rank 1 matrix, the conclusion is a direct consequence of 
theorem 2. Conversely if we know that V(g, Q') tr{MQMQ') = tr{MQ)tr{MQ'), 
then the linear form Q' 1-^ tr{MQMQ') — tr{MQ)tr{MQ') vanishes, which implies that 
MQM — tr{MQ)M for every matrix Q. Then, the theorem 2 tells us that M is a rank 1 
matrix. □ 
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